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ABSTRACT 

This paper reports positive detections of surface differential rotation on two rapidly 
rotating cool stars at several epochs, by using stellar surface features (both cool spots 
and magnetic regions) as tracers of the large scale latitudinal shear that distorts the 
convective envelope in this type of stars. We also report definite evidence that this 
differential rotation is different when estimated from cool spots or magnetic regions, 
and that it undergoes temporal fluctuations of potentially large amplitude on a time 
scale of a few years. 

We consider these results as further evidence that the dynamo processes operat- 
ing in these stars are distributed throughout the convective zone rather than being 
confined at its base as in the Sun. By comparing our observations with two very sim- 
ple models of the differential rotation within the convective zone, we obtain evidence 
that the internal rotation velocity field of the stars we investigated is not like that of 
the Sun, and may resemble that we expect for rapid rotators. We speculate that the 
changes in differential rotation result from the dynamo processes (and from the under- 
lying magnetic cycle) that periodically converts magnetic energy into kinetic energy 
and vice versa. 

We emphasise that the technique outlined in this paper corresponds to the first 
practical method for investigating the large scale rotation velocity field within convec- 
tive zones of cool active stars, and offers several advantages over asteroseismology for 
this particular purpose and this specific stellar class. 

Key words: Line: profiles - Stars: activity - Stars: imaging ~ Stars: binaries - Stars: 
magnetic fields - Stars: rotation. 



1 INTRODUCTION 

Despite considerable progress obtained over the last century, 
astrophysicists are still far away from a thorough under- 
standing on how magnetic fields are produced in the Sun and 
other cool active stars; the picture gets even worse when it 
comes to explain quantitatively how these fields generate the 
plethora of active phenomena observed both at their surface 
(e.g. dark spots) and in their immediate surrounding (e.g. 
corona, wind). Particularly interesting in this respect is the 
case of fully convective stars for which conventional dynamo 
processes (supposed to concentrate mostly in the interface 
layer between the radiative interior and the convective en- 
velope, i.e. the only place where the field can be stored for a 
time comparable to the period of the activity cycle) cannot 
be invoked to justify, not only the presence of a large-scale 
magnetic field (e.g. Johns-KruU & Valenti 1996), but also 
their very intense activity level. 



The most recent observational results in this field, con- 
sisting of temporal series of magnetic maps obtained by in- 
direct tomographic imaging from sets of rotationally mod- 
ulated spectropolarimetric data (Donati et al. 1992; Donati 
& Cameron 1997; Donati 1999; Donati et al. 1999; Donati et 
al. 2003; Petit et al. 2003a, 2003b), are particularly intrigu- 
ing and thus potentially very fruitful for providing new clues 
to what may be happening in the convective layers of these 
stars. All reconstructed magnetic maps to date indeed show 
the presence of large magnetic regions in the stellar photo- 
sphere where the field is mostly azimuthal (i.e. parallel to 
the surface and oriented along parallels); moreover, these re- 
gions are often observed to form partial, or even complete, 
rings of azimuthal field encircling the star at various lati- 
tudes, and are thus interpreted as the direct detection, at 
photospheric level, of the quasi-axisymmetric toroidal com- 
ponent of the large-scale dynamo field. 

These results suggest in particular that dynamos oper- 
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ating in very active stars are likely to be significantly differ- 
ent than that of the Sun. The observations indeed reveal that 
such dynamos are able to produce large scale field (and espe- 
cially toroidal field) very close to the surface (since we would 
not observe toroidal field in the photosphere otherwise) , and 
presumably even throughout the whole convective zone; dy- 
namo processes in these stars are therefore likely distributed 
within the convective envelope, rather than being confined 
at its base as in the Sun. This conclusion raises however 
a number of subsequent, and yet unresolved, problems; we 
can no longer understand, for instance, how the magnetic 
field can be stored within the convective zone for periods of 
time as long as decades. At the same time, it may provide 
new insight for a number of long standing mysteries, and in 
particular for our comprehension on how dynamos operate 
in fully convective stars. 

One of these mysteries concerns the long-term fluctua- 
tions observed in the orbital period of close binary stars in 
which one member is at least a cool active star (e.g. Hall 
1990), and occurring on a time scale of a few decades, i.e. 
much shorter than what would be required for tidal cou- 
pling between system components to operate (Zahn 1989). 
The least unreasonable interpretation for this phenomenon 
remains that proposed by Applegate (1992); he suggests that 
dynamo processes at work in one star of the binary system 
activate a periodic exchange between magnetic and kinetic 
energy within the convective zone of this star, and thus a 
cyclic fluctuation of its quadrupolar moment as well as of 
its gravitational fleld. Given the amplitude of the observed 
fluctuations in orbital period, one can easily infer that such 
magnetic to kinetic energy exchanges speciflc to dynamo 
processes must occur throughout the whole convective zone 
and not only within a thin layer (Lanza et al. 1998; Do- 
nati 1999), i.e. that dynamo processes in these stars should 
indeed be distributed. 

If Applegate's (1992) idea is true, it therefore implies 
that very active stars such as those found in close binary 
stars, but also young single stars that did not dissipate most 
of their angular momentum yet and still exhibit very ener- 
getic active phenomena, should feature a convective zone 
that globally undergoes such periodic exchanges between 
magnetic and kinetic energy. It means in particular that 
the internal velocity field within the convective zone of very 
active stars, and therefore both their radial and surface dif- 
ferential rotation profiles, should vary with time in a way 
correlated to the magnetic cycle. Detecting such variations 
in the differential rotation of active stars thus appears as a 
very interesting observational challenge, as it would provide 
a definite confirmation of both Applegate's (1992) mecha- 
nism and of the existence of distributed dynamos. It would 
also bring into a much wider astrophysical context the re- 
cent discovery that the solar angular rotation at the base 
of the convective zone is undergoing temporal fluctuations 
(Howe et al. 2000) as a probable consequence of the activity 
cycle. 

Measuring surface differential rotation of very active 
stars is however a rather tricky task; detecting changes of 
this differential rotation is even more difficult. Several meth- 
ods have been proposed to estimate surface differential ro- 
tation. Some of them try to make use of the subtle changes 
that differential rotation induce in the profile of spectral 
lines (Brunning 1981; Reiners & Schmitt 2002); this method 



can however not be used for the stars we are interested in, 
the spectral lines of which being heavily distorted by the 
presence of cool spots at their surfaces. Other methods (e.g. 
Donati & Cameron 1997) propose to take advantage of these 
surface features, and use them as tracers to derive some in- 
formation on how their rotation periods depend on latitude. 
The most recent of such techniques, that of Cameron et al. 
(2002), succeeded, not only in estimating surface differen- 
tial rotation of one active star, but also in suggesting that 
temporal changes in the amount of differential rotation were 
indeed occurring in this star (Cameron & Donati 2002). 

In this paper, we propose to make use of yet another 
technique, first proposed by Donati et al. (2000) to estimate 
the differential rotation of a young pre-main-sequence star, 
then by Petit et al. (2002, 2003a, b) both for simulation 
purposes and application to spectropolarimetric data. The 
aim of this paper is to apply this method to the extensive 
spectropolarimetric data set that we collected in the last 
seven yr for three stars, namely the young ultra-fast rotator 
AB Doradus, the young KO star LQ Hydrae and the Kl sub- 
giant of the RS CVn system HR 1099, and from which yearly 
brightness and magnetic maps were obtained and published 
in the literature (Donati & Cameron 1997; Donati 1999; Do- 
nati et al. 1999, 2003). In Sect.|21 we recall the main aspects 
of both observational material and modelling tool, then de- 
scribe its application to the three selected stars in Sect.|21 
after discussing at length the implication of our results for 
the understanding of the global dynamics of convective zones 
of cool stars (in Sect.|lJ, we finally conclude and propose in 
Sect.|S]a few directions in which this work could be fruitfully 
extended. 



2 OBSERVATIONAL MATERIAL AND 
MODELLING TECHNIQUES 

The data we use are the spectropolarimetric observations 
collected at the Anglo-Australian Telescope (AAT) with a 
visitor polarimeter mounted at Cassegrain focus and fibre 
linked to the high resolution UCL Echelle Spectrograph 
(UCLES). This material, and in particular the observing 
logs, the observing procedures and the data reduction de- 
tails, are described extensively in a series of published pa- 
pers that present the results obtained up to now (Donati et 
al. 1997; Donati & Cameron 1997; Donati 1999; Donati et 
al. 1999, 2003). 

The data we use in this paper are all observations of 
AB Dor, LQ Hya and HR 1099 collected since epoch 1995.9, 
i.e. obtained in a similar spectrograph configuration. This 
makes a maximum amount of 7 data sets (corresponding to 
epochs 1995.94, 1996.99, 1998.03, 1999.00, 1999.97, 2000.93 
and 2001.99) for each of the three stars; we however antici- 
pate that several of them (and in particular those with very 
few data or collected over a timespan of no more than a few 
days, e.g. those at epoch 1999.97) will be unusable for our 
purpose, according to the conclusions of Petit et al. (2002) . 
Note that all data were reprocessed with the newest version 
of the reduction software (Donati et al. 2003) to ensure the 
highest possible internal accuracy in radial velocity as well 
as the best level of homogeneity. All data were run through 
Least-Squares Deconvolution (LSD), a multiline tool devel- 
oped by Donati et al. (1997) found to extract successfully 
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Figure 1. Reduced maps in the differential rotation parameter plane, obtained for AB Dor from the Stokes I (left panel) and Stokes 
V (right panel) data sets collected at epoch 2001.99. Black corresponds to models for which the minimum level was obtained (given 
the assumed reconstructed information content), while white depicts those for which the is larger by at least 1% (i.e. by least 274 and 
69 for the left and right panels respectively) than the minimum x^- The outer envelope of all grey points corresponds to 11 and 5.5 a 
confidence ellipses for both differential rotation parameters (for the left and right panels respectively), while its projection on each axis 
gives the 16 and 8 cr confidence intervals for each parameter taken separately. Note that distances on both axes are not equal. 



the average rotationally broadened shape of unpolarised or 
circularly polarised lines (called Stokes / and V LSD pro- 
files in the following) at each observed phase of the stellar 
rotational cycle, from the thousands of moderate to strong 
spectral features present in the recorded wavelength domain. 

As mentioned earlier, the method we chose to estimate 
surface differential rotation for the selected stars is that in- 
troduced by Donati et al. (2000), the validity of which was 
confirmed through an extensive sets of simulations by Pe- 
tit et al. (2002). Recently, Petit et al. (2003a, b) applied it 
successfully to large data sets collected over several years on 
a different telescope and with a different (though similar) 
instrumentation for two different stars, one of them being 
our third object of interest (the Kl subgiant of HR 1099) 
and the other a giant star of FK Com type (HD 199178). 
This method consists in assuming a given surface differen- 
tial rotation law (that we implement in our indirect stellar 
surface imaging code) and producing brightness or magnetic 
images of the stellar surface from all unpolarised or circu- 
larly polarised data available for this star at a given season; 
we then take as most probable the differential rotation law 
that generates the images with lowest information content 
(given a certain data fit accuracy), or equivalently the im- 
ages that provide the best fit to the data (given a certain 
image information content). In practice, both options yield 
very similar results; the second one happens to be more con- 
venient as it provides an easy way of estimating error bars on 
the differential rotation parameters (see Petit et al. 2002). 

The stellar surface imaging code that we use is that of 
Brown et al. (1991) and Donati & Brown (1997), applied 
quite extensively to real data in the last decade (the latest 
examples being Donati et al. 2003 and Petit et al. 2003a, 
b). It includes in particular the possibility of reconstructing 
images of stars that are differentially rotating, i.e. of taking 
into account explicitly the differential rotation pattern in 
the imaging process. The surface differential rotation law 
we assume in this paper is given by the following equation: 



Q{e) = 0,^^- dfl cos^ e (1) 

where Q{8) is the rotation rate at colatitude 9, Qcq the rota- 
tion rate of the equator and dQ the difference in rotation rate 
between the pole and equator. At each season, one bright- 
ness and one magnetic image is produced from the full Stokes 
I and Stokes V data sets respectively, for each selected pair 
of differential rotation parameters within a reasonable search 
domain of the Qcq-di} plane. This typically means that we 
have to compute a few hundred brightness and magnetic im- 
ages for each star at each season, which therefore requires 
large amount of computing time; for instance, in the par- 
ticular case of the Stokes I data set of AB Dor at epoch 
2001.99 where as much as 27,400 data points must be fitted 
simultaneously, we typically need a few days of computing 
time on the fastest available workstations to derive how the 
reduced of the fit to the data varies with the differential 
rotation parameters at a good enough resolution. A typical 
example of the reduced maps we obtain from Stokes / 
and V data is shown in Fig. 

It could be argued that, in the sake of computing time, 
we could actually derive the optimal differential rotation 
parameters as part of the imaging process itself, since this 
would mean no more than fitting another two parameters in 
addition to the thousands of image pixels we already have to 
deal with. However, our experience is that the maps we 
obtain are not always as nice as those shown in Fig. Q and 
sometimes do not feature any physically meaningful mini- 
mum. We therefore preferred to keep our initial procedure, 
and check in each case whether the obtained reduced x^ 
maps showed a clear enough minimum to be able to claim 
safely that differential rotation is indeed detected at the 
corresponding epoch. Another very important advantage of 
producing reduced x^ maps is that they can also be used 
in a very straightforward way to obtain error bars on our 
estimates of the differential rotation parameters. Following 
Press et al. (1992), we indeed know that the 1 a confidence 
interval for each differential rotation parameter taken sep- 
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arately can be obtained by searching for all points of the 
f2eq-dn plane for which the increase (with respect to the 
minimum value in the map) is smaller than 1^ (or equiv- 
alently for which the reduced increase is smaller than 
1/n, where n is the number of fitted data points), and pro- 
jecting this ellipsoid onto the f2eq and dfi axes respectively. 
The simplest way to evaluate these projections (and thus 
the formal error bars for each parameter) is to compute the 
curvature radii of the paraboloid at its minimum and 
the correlation coefficient between the two parameters; for 
a paraboloid with linear and quadratic terms given by 
a Slcq^ + h fioq dO. + c + d ficq + e dfi, the value of 
the optimal parameters and the associated error bars sim- 
ply write: 

^cq = (te-2cd)/(4ac-6^), cT{Q.r^) = [a - / -icY^''^ (2) 

dn = {hd~2ae)/{Aac~b^), cj{dQ) = [c - / AaY^''^ (3) 

The procedure we thus carry out to obtain both the op- 
timal differential rotation parameters and associated error 
bars simply consists in fitting the reduced maps (actually 
only the points in the neighbourhood of the minimum) by a 
bi-dimensional paraboloid (as expressed above), derive the 
five coefficients of this paraboloid and then the requested 
quantities using Eqs. Inland |3] In addition to being the eas- 
iest way to evaluate both the position of the minimum and 
the curvature radii at that minimum, it also averages out 
the slight random differences in convergence accuracy that 
always remain at a relative level of 0.01% between all points 
in the map. We find this method to be very robust in prac- 
tice, and the derived differential rotation parameters very 
weakly sensitive to internal parameters, such as for instance 
the amount of image information towards which the code is 
asked to converge. 



3 SURFACE DIFFERENTIAL ROTATION 

As explained in Petit et al. (2002), estimating the surface 
differential rotation of a cool spotted star requires that the 
spectroscopic (or spectropolarimetric) data set used in this 
aim correctly samples in time the variability in the sur- 
face distribution of cool spots or magnetic regions associ- 
ated with this phenomenon. It means in particular that 
large fractions of the stellar surface must be observed several 
times, with adequate temporal gaps in between. It straight- 
forwardly tells us that sparse data sets, such as those ob- 
tained for all three stars in 1998 January (epoch 1998.03) 
and 1999 December (epoch 1999.97) in which most regions 
of the stellar surfaces was observed only once, are not be 
appropriate for our purposes. Similarly, the data set col- 
lected on LQ Hya at epoch 1995.94, containing only two 
main groups of observations (plus a third one of very poor 
quality, Donati et al. 1997; Donati 1999), is totally useless 
for the present investigation. 



These error bars should not be confused with the error ellipses 
for both parameters taken simultaneously, such as those depicted 
in Figs.|21to|l] which includes all points of the map for which the 
X'^ increase is smaller than 2.3 



3.1 The young ultra-fast rotator AB Doradus 

As far as AB Dor is concerned, it leaves us with five 
main data sets (corresponding to epochs 1995.94, 1996.99, 
1999.00, 2000.93 and 2001.99) with adequate phase over- 
lap between different groups of observations. The main rea- 
son for this is of course that the rotation period of AB Dor 
is short enough (0.51479 d, equivalent to a rotation rate 
of 12.2053 radd~^) to allow a large fraction of the stellar 
surface (as much as two thirds) to be observed in a single 
night, and sufficiently close to half a day that reasonable 
phase overlap is obtained between data sets separated by a 
few days. The estimates we obtain are listed in Tableland 
shown visually on Fig. |5| The imaging parameters that we 
assumed for AB Dor in this experiment are those derived by 
Donati et al. (2003), with the inclination of the rotation axis 
to the line of sight i set to 60° and the projected equatorial 
rotation velocity vsaii set to 89 kms~^. 

Several comments can be made from these results. First 
of all, we want to stress the overall good agreement of the 
estimates obtained from either the Stokes I or the Stokes 
V data. If the subtle signature we are attempting to re- 
trieve from our spectra was not produced by surface dif- 
ferential rotation, but rather by some spurious problem in 
the analysis, then there would be no reason for us to find 
roughly the same result five times, from five totally inde- 
pendent data sets. It represents already a strong argument 
that what we see is indeed surface differential rotation. We 
can then also note that the values we derive from Stokes 
/ data at epoch 1995.94 are in good agreement (i.e. within 
about one error bar) with the estimates we originally pub- 
lished (r2cq = 12.2434 radd"\ dn = 56.4 mradd"^) when 
differential rotation was first detected on AB Dor (Donati 
& Cameron 1997), indicating that the new method we are 
using is well behaved and provides results that are compat- 
ible with those yielded by older and more straightforward 
(though less accurate and rigorous) techniques. 

There is however a number of significant differences be- 
tween the results obtained recently with yet another tech- 
nique (estimating the recurrence rate and latitude of indi- 
vidual starspots by accurately tracing the repetitive trails 
they produce in dynamic spectra each time they cross the 
visible stellar hemisphere) and published by Cameron & 
Donati (2002). For the two epochs (actually the two data 
sets) we have in common, these authors obtain differen- 
tial rotation parameters of Qcq ~ 12.2498 ± 0.0019 radd^^ 
and dn = 71.4 ± 5.8 mradd"^ for epoch 1995.94, and 
ficq = 12.2417±0.0017radd"^ and dQ = 57.6±4.8 mrad d'^ 
for epoch 1996.99, i.e. way out both our and their sets of 
errors bars. Note that the variation of both parameters be- 
tween the first and second epochs is found to be roughly 
similar in both studies, suggesting that most of the discrep- 
ancy is likely attributable to slight internal calibration errors 
in the second method but does not question the reality of 
the measurements themselves. We indeed find that this sec- 
ond method is much more sensitive than the one used here to 
the assumed v sin i and that part (about half) of the discrep- 
ancy we report comes from the different v sin i used in both 
studies, set to 89 kms~^ for our experiment (as obtained by 
Donati et al. 2003) and to 91 kms~^ in Cameron & Donati 
(2002). The remaining difference seems to be attributable to 
the way Cameron & Donati (2002) weigh the rotation rates 
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Table 1. Summary of differential rotation parameters obtained for AB Dor on each season. Columns 2 to 6 correspond to Stokes / data 
while columns 7 to 11 to Stokes V data. Cols. 2 and 7 list the derived equatorial rotation rate Oeq with its 68% (i.e. 1 cr) confidence 
interval, Cols. 3 and 8 the difference in rotation rate dQ between the equator and pole with its 68% confidence interval, Cols. 4 and 9 
gives the inverse slope of the ellipsoid in the Ocq— rfH plane (also equal to cos'^ 9s where 9^ denotes the colatitude of the gravity centre of 
the spot distribution, see Donati et al. 2000), Cols. 5 and 10 the rotation rate Qs at colatitude 9s, and Cols. 6 and 11 the total number 
of data points used in the imaging process. 
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Figure 2. Differential rotation parameters of AB Dor as measured from Stokes / (left panel) and Stokes V data (right panel). All points 
are marked with + signs and joined together with a dashed line by order of increasing epoch starting from the first point in the series 
(corresponding to epoch 1995.94 and marked as o). For each point, 68% confidence ellipses for both parameters (projecting on both axes 
into 1.5 a confidence intervals for each parameter taken separately) are also depicted. Note the different scale used for each plot. 



derived from individual spots when it comes to fitting them 
with a differential rotation law; if, instead of using roughly 
equal weight for all spots as they did, larger weight is at- 
tributed to bigger spots (as our technique implicitly does), 
a much better (though still not quite perfect) agreement is 
obtained. It is not clear yet whether this problem reflects no 
more than numerical uncertainties, or something intrinsic to 
the star, e.g. that small and large spots do not happen to 
suffer the same latitudinal shear (as is the case on the Sun) . 

Although the different results obtained from either 
Stokes / and V data are roughly all of the same magnitude, 
they nevertheless show definite differences among them, as 
one can see quite straightforwardly from Fig.|5| From Stokes 



I data for instance, we note that f^oq changed by as much 
as 6.2 mradd"^ between epoch 1995.94 and 1996.99, about 
7.8 times the quoted error bar; comparing now how dO, 
evolved between epoch 1996.99 and 1999.00, we obtain a 
variation of 11.2 mradd"^, about 4.5 times the largest error 
bar. Similarly, results from Stokes V data indicate variations 
as large as 13.7 and 37.1 mradd"^ in f2eq and dTl respec- 
tively, again about 4.5 times the quoted error bars (of 3.1 and 
8.1 mradd"'^ respectively). This result by itself is enough to 
safely claim that we indeed detected fluctuations in the dif- 
ferential rotation parameters of AB Dor (with a false alarm 
probability much smaller than 1 ppm), thus confirming the 
initial report of Cameron & Donati (2002). Our conclusion 
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relies of course on the fact that our error bars are not sig- 
nificantly underestimated; although the simulations of Petit 
et al. (2002) reveal that this may indeed happen when the 
phase sampling of the observations is sparse, they neverthe- 
less establish quite firmly that this does not occur with as 
dense a sampling as that secured for AB Dor. 

A rather surprising discovery is that the differential ro- 
tation parameters obtained from Stokes V data are clearly 
larger, by an amount that varies between 25% to as much 
as 80%, than those derived from Stokes / data, for a reason 
that we do not fully understand yet. The effect in itself is 
already clearly visible on the raw maps themselves, such 
as those shown in Fig.Qfor instance. One can in particular 
wonder whether it is not due to some artifact inherent to the 
method itself. The first natural suspicion is that the surface 
differential rotation law we assumed in the imaging process 
is significantly different from the true one; in this case, po- 
tential differences in the latitude range respectively sampled 
by the brightness and magnetic regions recovered at a given 
epoch may explain the observed discrepancy between values 
of r2cq and dO. derived from Stokes I and V data. If this was 
the correct interpretation however, we would then expect to 
see larger discrepancies when the latitudinal barycentres of 
the brightness and magnetic distributions are further away 
from one another. This is certainly not what we get, e.g. at 
epoch 2000.93 where both the latitudinal barycentres of the 
brightness and magnetic distributions almost coincide (see 
Table while the photospheric shear derived from Stokes 
V data exceeds that deduced from Stokes / data by more 
than 60%. We thus have to conclude that the origin of the 
observed discrepancy is somewhere else. 

Another potential explanation may be looked for in the 
fact that magnetic features hosting radial field lines, suf- 
fering from a stronger dependence to limb darkening than 
brightness features (Brown et al. 1991), are reconstructed 
with a stronger bias towards higher latitudes than cool spots; 
matching the observed recurrence rate of surface features 
would then require larger S7cq and dSl when magnetic re- 
gions are considered rather than cool spots. This however 
cannot be the case. The differential latitude bias that we 
mention remains indeed rather small, especially when phase 
coverage is as dense as that we have on most seasons, and 
does certainly not exceed a few degrees. Moreover, obser- 
vations indicate that magnetic field is, most of the time, 
predominantly azimuthal (Donati et al. 2003) rather than 
radial, and azimuthal field regions are actually not sensitive 
to the above mentioned effect. We can thus safely claim that 
the considerably larger shear observed from Stokes V data 
cannot be explained this way. Again, the best method for 
investigating whether this phenomenon is genuine or spu- 
rious is to carry out simulations; this is what Petit et al. 
(2002) did, leading to the conclusion that no such bias was 
to be observed between differential rotation parameters es- 
timated from either Stokes / or Stokes V data when phase 
sampling was as dense as ours. Similar simulations, either 
on fake data or on the real observations presented here, also 
enable us to conclude that reasonable errors in the imaging 
parameters (for instance errors of the order of 2 kms~^ in 
II sin i) produce only small variations in the results (of the 
order of one error bar typically) and in particular negligible 
changes in this discrepancy we report between differential 
rotation parameters obtained from Stokes / and V data. 



The only conclusion we are left with is that the ob- 
served difference reflects some genuine property of surface 
features in cool active stars, for instance that magnetic re- 
gions are not anchored at the same depth in the convective 
zone and thus do not witness the same shear. Direct evidence 
from this can be found by comparing the average latitude at 
which cool spots and magnetic regions are found to appear, 
as well the estimated recurrence rate at this average latitude 
(both parameters being listed in Table^. At epoch 1999.00 
for instance, we obtain that magnetic regions preferentially 
form at a latitude of 39.9° and rotate at an average rate of 
12.219 radd~^, while cool spots tend to show up at a lati- 
tude 23.3° and rotate at a mean rate of 12.230 radd~^; the 
latitudinal shear we derive from this, equal to 43 mradd"^, 
is significantly smaller than that obtained by using either 
Stokes / data alone (58.4 ± 1.5 mradd"^) or Stokes V data 
only (79.5 ± 7.4 mradd"^). It clearly indicates that both 
sets of features do not refer to the same latitudinal angu- 
lar velocity field, and therefore likely not to the same depth 
in the convective zone. We will come back on this point in 
Sect.0] Note as well that the pattern of temporal variation 
in the differential rotation parameters is not the same either 
between Stokes I and Stokes V data; looking at each change 
from one season to the next in the Q.cq-d^ plane as a vector 
of given length and direction, we observe that the vectors 
associated with Stokes I data do not necessarily feature the 
same direction as (and sometimes even point in the opposite 
direction to) those associated with Stokes V data (the most 
obvious example being the evolution between epoch 2000.93 
and 2001.99). Again, we consider this as additional evidence 
that both sets of features correspond to different zones of 
the convective envelope. 

Note that this difference in differential rotation as re- 
vealed by cool spots and magnetic features was not seen 
in the initial study of Donati & Cameron (1997), probably 
because only one component of the reconstructed magnetic 
field (the radial field map) was used at that time in the 
cross correlation experiment from which the photospheric 
shear was measured; in the present study, we make use of 
all available information in Stokes V profiles, implying that 
the present results are definitely more reliable. 

3.2 The young KO dwarf LQ Hydrae 

For LQ Hya, only three data sets (those corresponding to 
epochs 1996.99, 2000.93 and 2001.99) happened to be rich 
enough to yield a positive detection of surface differential ro- 
tation. The one obtained at epoch 1999.0, although featuring 
a rather large number of observations (54 profiles in Stokes 
/ and 14 in Stokes V) yields maps with no clearly defined 
paraboloid, most likely as a result of the much poorer qual- 
ity of the data secured during this campaign (see Donati et 
al. 2003). The opposite situation is observed with the data 
collected at epoch 1996.99; although gathered in a period 
of only 5 d, they are nevertheless able to provide us with a 
meaningful estimate of surface differential rotation, but only 
from the Stokes V subset. The results we obtain are listed 
in Table |5| and shown in Fig.|21 The imaging parameters we 
use for LQ Hya are again those determined by Donati et al. 
(2003), equal to i = 60° and vsm.i = 26 kms~^. 

The first thing worth reporting is of course the detec- 
tion of differential rotation itself at the surface of LQ Hya; al- 
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Table 2. Same as Table for LQ Hya 







Stokes / data 








Stokes V data 






Epoch 






cos^ 9s 
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cos^ 9s 
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(mrad d^-"^) 


(mradd^-'^) 




(radd-1) 




(mradd"-'^) 


(mradd"^) 




(radd-i) 




1996.99 












3, 892.3 ±9.6 


-48.7 ±23.5 


0.355 


3.910 


299 


2000.93 


3, 982.4 ±9.1 


194.2 ± 21.6 


0.398 


3.905 


2,835 


3, 969.3 ±22.9 


200.9 ±72.4 


0.298 


3.909 


805 


2001.99 


3, 922.3 ±1.0 


14.4 ± 2.9 


0.333 


3.918 


5,945 


3, 944.6 ±2.3 


86.7 ±8.7 


0.213 


3.926 


1,558 
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Figure 3. DifTerontial rotation parameters of LQ Hya. The 
dashed line links together measurements from Stokes V data by 
order of increasing epoch (starting from epoch 1996.99). As in 
Fig-IH 68% confidence ellipses for both parameters are also shown. 



though this discovery is not altogether very surprising (since 
we probably expect all stars with a significant convective 
envelope to rotate differentially), it is nevertheless the first 
detection of differential rotation on this star and the second 
one on a zero-age- main-sequence star (after AB Dor), mak- 
ing it in this respect very interesting for potential compar- 
isons with the results obtained on AB Dor. In that respect, 
the very large photospheric shear detected at epoch 2000.93 
(equal to about 0.2 radd~^, implying a time for the equator 
to lap the pole by one complete turn of only one month) is 
about 3±1 times larger than that of the Sun or than the aver- 
age dO, we measure for AB Dor. Given the fact that LQ Hya 
and AB Dor are very similar in terms of internal structure 
(Donati et al. 2003), it already clearly demonstrates that 
any relation we may eventually derive from the observations 
and giving the strength of this shear as a function of, say, 
stellar mass (e.g. Kitchatinov & Riidiger 1999) is expected 
to show a significant scatter. 

One possible origin for this scatter is of course potential 



variability with time, and this is precisely what we observe 
for LQ Hya, with a shear varying from 0.20 radd~^ at epoch 
2000.93 to only 0.014 and 0.087 radd"^ for Stokes I and V 
data respectively at epoch 2001.99, i.e. within only one year. 
Note that this variation of dfi as measured from Stokes / 
data is detected at a level of 8.3 a while the corresponding 
one in f7cq amounts to about 6.6 a, implying that both can 
be considered as fully reliable. A similar, though inverse, 
variation must have occurred at least once between epoch 
1996.99, at which we detect no differential rotation (or at 
least a value compatible with solid body rotation within 2 a) 
at an accuracy of 0.024 radd~^, and epoch 2000.93 where 
this very strong shear of 0.2 radd~^ is observed. It suggests 
that the convective envelope may indeed be passing cycli- 
cally from stages where it rotates almost as a solid body, to 
others where it exhibits strong differential rotation between 
the equator and pole (and likely throughout its depth as 
well); we discuss this point more extensively in Sect. 3] 

As for AB Dor, we also notice a clear difference between 
the shear as estimated from Stokes / and V data. While the 
shear measured from unpolarised data at epoch 2001.99 is 
only 0.014 radd"\ the one derived from Zeeman signatures 
is about 6 times larger! Note that larger shears are again 
obtained from magnetic images, as was already the case for 
AB Dor. No such difference is observed at epoch 2000.93; one 
must however keep in mind that error bars at that epoch are 
considerably larger than those from our most recent data set. 
We again speculate that this difference most likely refiects 
the fact that the cool spots and magnetic regions we observe 
in our reconstructed images are dissimilar in their essence, 
and turn out in particular to be anchored at different depths 
in the convective zone. We must admit that this conclusion 
is very surprising, although not totally unexpected from the 
fact that only very weak spatial correlation was observed 
between the reconstructed magnetic and brightness features 
(e.g. Donati & Cameron 1997; Donati 1999; Donati et al. 
2003); the fact that these two classes of photospheric fea- 
tures do not witness the same rotational shear strengthens 
the reality of this difference, and may eventually give us 
more insight on its physical nature. 



3.3 The Kl subgiant of the RS CVn system 
HR 1099 

For HR 1099, only two data set turned out to be of interest 
for our particular purpose; while those obtained on 1995.94 
and 1996.99 are essentially too sparse to yield anything use- 
ful, that secured at epoch 1999.0 is plagued by numerous 
observations with poor quality (see Donati et al. 2003), as 
was already the case for LQ Hya. Even Stokes I data from 
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Table 3. Same as TableE]for HR 1099 



Stokes / data Stokes V data 

Epoch r2eq dQ. cos^ 9s Sis n Slcq do, cos^ 9s f!s n 

(mradd"-'^) (mradd^-'^) (radd^-"^) (mradd~^) (mradd^^) (radd^-"^) 

2000.93 2, 222.2 ±4.2 21.2 ±8.1 0.477 2.212 1,519 

2001.99 2, 223.5 ±2.1 -2.5 ± 4.8 0.425 2.225 6,156 2, 220.3 ±3.2 -0.4 ±7.1 0.431 2.221 1,539 



HR 1099, Stokes I &l V data 

n ' \ ' I 




2.2 2.22 2.24 

Oeq (rad/d) 

Figure 4. Same as Fig. Elfor HR 1099. Tlie point marked witii 
o now corresponds to epoch 2000.93. 



epoch 2000.93 are found to produce maps with no clearly 
defined minimum, but rather an infinitely long valley char- 
acterising cases where our method fails. The exact reason 
for which this failure occurs is not entirely clear yet; all con- 
ditions for a success of the method, and in particular the 
presence of cool spots on the brightness image at both low 
and high latitudes as well as a relatively long series of obser- 
vations that regularly come back on the same regions of the 
stellar surface, were however met in this data set. We specu- 
late that it may be due to the fact that some other source of 
variability in the spot distribution, intrinsic in nature (like 
for instance the formation of new spots and the disappear- 
ance of old ones), was operating at a much higher level at 
that specific epoch, preventing the signal from the system- 
atic evolution induced by differential rotation to build up 
properly in the maps. The results we obtain are listed 
in Table |3 and shown in Fig. |1] and were derived using 
the imaging parameters determined by Donati et al. (2003) , 
equal to i = 38° and usini = 39 kms"'^. 

The error bars we obtain are still too large to claim 
that HR 1099 is not rotating as a solid body; even at epoch 
2000.93 where the signal is strongest, the value of dfi we de- 



rive differs from by no more than 2.6 a, which we do not 
consider as large enough to be called a reliable detection. 
Our estimates nevertheless indicate altogether a tendency 
for the differential rotation of HR 1099 to be solar like (i.e. 
with a positive dO.), at least at the time of our observations, 
rather than antisolar (as announced by various authors in 
the recent years, e.g. Vogt et al. 1999; Strassemeier & Bartus 
2000). Our measurements also pin down the strength of the 
photospheric shear to something of the order of 10 mrad d~^, 
in reasonable agreement to the actual estimates obtained by 
Petit et al. (2003a) from independent data sets. If we aver- 
age together the values of dO, derived from Stokes / and V 
data at epoch 2001.99, we obtain that the shear is compati- 
ble with 0, with an error bar of 5 mradd"^. This is slightly 
smaller than the value found by Petit et al. (2003a) at the 
same epoch (equal to 18±4 mradd"^) and may indicate that 
our error bars are slightly underestimated (as suggested by 
the simulations of Petit et al. 2002 in a case with similar 
phase sampling). Since no clear detection is obtained, noth- 
ing can of course be quoted on a potential variation with 
time of the differential rotation parameters. 



4 SUMMARY AND DISCUSSION 

4.1 Differential rotation, dynamo processes and 
activity cycles 

We applied in this paper a new technique for measuring sur- 
face differential rotation in cool active stars featuring bright- 
ness and magnetic spots on their surface, and using these 
surface features as tracers of the stellar rotation at various 
latitudes. First introduced by Donati et al. (2000) then val- 
idated through extensive simulations by Petit et al. (2002), 
this method is found to produce very accurate and reliable 
results on the three stars we applied it to in this new study. 

We obtain positive detections of surface differential ro- 
tation in two of our three stars (namely the young KO dwarfs 
AB Dor and LQ Hya). In the first of these two stars, the pho- 
tospheric shear (and equatorial rotation rate) we measure is 
in very good agreement with earlier estimates published by 
Donati & Cameron (1997) and Donati et al. (1999) from the 
same data sets, but with a different analysis. In the second 
object, surface differential rotation is detected for the first 
time, and turns out to be rather large (as much as 3±1 times 
that of the Sun at epoch 2000.93). In both stars and at all 
epochs at which our method is successful, the photospheric 
shear that we measure is found to be solar like, i.e. with the 
stellar equator rotating faster than the pole. No evidence 
for anti solar differential rotation such as that claimed by, 
e.g., Vogt et al. (1999) or Strassemeier & Bartus (2000) is 
found in our data. For the third star (the Kl subgiant of 
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HR 1099), we only obtain that the photospheric shear (il 
present) is of the order of 10 mradd"^ (i.e. about 5 times 
weaker than that of the Sun), in good agreement with the 
recent estimates obtained by Petit et al. (2003a). 

For both stars on which differential rotation is detected, 
we also observe significant variation of the differential rota- 
tion parameters with time. This confirms in particular the 
earlier claim by Cameron & Donati (2002) that the pho- 
tospheric shear of AB Dor was fluctuating from epoch to 
epoch, and amplifies the importance and significance of this 
discovery by revealing that other similar stars actually be- 
have in a way roughly identical in essence, though much 
more extreme in intensity. Our observations indeed disclose 
that, while the temporal variation of the differential rota- 
tion parameters remains moderate in AB Dor (smaller than 
40 mradd"^ for dfl and 15 mradd"^ in Q^q peak-to-peak 
throughout our 5 epochs of observations), it can be much 
larger in LQ Hya where it can switch from virtually solid 
body rotation (e.g. at epoch 1996.99) to a shear as strong 
as 200 mradd"^ with an increase of about 80 mradd"^ in 
fieq, an effect about 5 times larger in amplitude than what 
we witnessed for AB Dor. 

Another important discovery reported in this paper is 
the fact that the differential rotation parameters we obtain 
from Stokes / and Stokes V data are not compatible to 
within the measurement error bars. We find that this effect, 
although strange, cannot be attributed to a systematic prob- 
lem in the method (see in particular Sect. 13. ll where this is 
discussed extensively). We are thus left with the conclusion 
that this discrepancy actually reflects some real property of 
brightness and magnetic surface features of active stars, and 
therefore of the convective zones in which they are anchored. 
A possible interpretation of this result (already mentioned 
above) is that the discrepant angular velocity laws respec- 
tively derived from magnetic and brightness features corre- 
spond to different layers within the convective zone, e.g. to 
different anchoring depths of brightness and magnetic sur- 
face regions. If this is the true explanation, the large dif- 
ference in angular shear witnessed by both types of surface 
features at a given epoch (e.g. reaching up to 72 mradd"^ 
between the pole and equator on LQ Hya at epoch 2001.99) 
strongly argues in favour of dynamo processes, and thus 
potential anchoring locations of surface features, being dis- 
tributed throughout the whole convective zone rather than 
being confined at its base as in the Sun. In this context, the 
observed discrepancy may serve as an indirect tool for study- 
ing angular velocity fields within stellar convective zones (see 
Sect. 

The first suspicion is that the temporal variations in 
differential rotation that we detect only reflect the fact that 
differential rotation is not a strict rule that all surface fea- 
tures must obey, but only an average behaviour that cool 
spots and magnetic regions follow in a statistical way; this 
is in particular what is observed on the Sun (with various 
classes of differential rotation tracers behaving differently 
and various tracers within each class departing from the av- 
erage differential rotation law for this class), reflecting the 
fact that differential rotation is a consequence of turbulence 
and inevitably contains some randomness when looking at 
individual features. However, the temporal variations we re- 
port here (especially in the case of LQ Hya) are so large that 
we simply cannot invoke the random nature of differential 



rotation to account for the observed variations, and have no 
better explanation to put forward than assuming that the 
average differential rotation law itself is changing as a whole 
from year to year. 

In this context, a natural idea is that these variations 
are caused by the stellar magnetic cycle converting period- 
ically kinetic energy within the convective zone into large 
scale magnetic flelds and vice versa, as originally proposed 
by Applegate (1992) in the particular case of close binary 
systems; if his scenario is correct, there is of course no rea- 
son for which all similarly active stars should not exhibit 
the same phenomenon. Although this interpretation would 
be conceptually very satisfying, it is nevertheless premature 
to conclude already that this is indeed what we observe; a 
definite demonstration would require for instance to mon- 
itor a few stars for a long enough time and see both the 
differential rotation parameters and the global polarity of 
the large scale dynamo field to show cyclical variations on 
the same temporal period, or at least to exhibit strongly 
correlated fluctuations in case of a non-cyclic behaviour. All 
activity proxies estimated from the reconstructed stellar im- 
ages (such as the total spot coverage or the mean quadratic 
magnetic flux) and monitored during the last decade did ex- 
hibit no clear evidence of cyclic variations with time (Donati 
et al. 2003a). 

4.2 Angular momentum and kinetic energy in a 
differentially rotating convective envelope 

Looking at how fieq varies as a function of dfl in any active 
star for which both parameters, as well as their variations, 
can be estimated can also provide us with additional clues on 
how angular rotation is distributed within their convective 
envelope. In this initial investigation, we suggest to explore 
this idea by assuming a priori speciflc angular velocity fields 
within the convective envelope, and derive what relation 
angular momentum conservation imposes between fieq and 
dn. One may argue that angular momentum conservation 
does not apply to the convective zone as a whole since this 
zone regularly exchanges angular momentum both with the 
circumstellar environment (through magnetospheric confine- 
ment of massive prominences and winds) and with the ra- 
diative interior (through a weak magnetic coupling for in- 
stance). These external phenomena are however known to 
occur on a much longer timescale than the redistribution 
of angular momentum within the convective zone itself; for 
the purpose of this first basic study, we can therefore safely 
assume that the convective zone obeys global conservation 
of angular momentum. 

The specific angular velocity fields we propose to con- 
sider in this paper are that of the Sun (with angular rotation 
independent on the radial distance to the centre of the star) , 
and that supposed to occur in very rapid rotators (with 
angular rotation constant along cylinders symmetric with 
respect to the rotation axis). Of course, we expect many 
more situations to occur potentially in real stars, but the 
two selected ones likely represent the extreme cases we may 
encounter among all active stars. In both cases, we assume 
that the surface differential rotation obeys Eq.0 We obtain 
that the internal angular velocity field n(r, 6) is given, in 
the first of our two test cases, by: 

n{r, 9) = fipoio + dD. sin^ 9 (4) 
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where r denotes the radial distance to the centre of the star, 
9 the colatitude with respect to the rotation axis and f2poic 
the rotation rate at the pole (equal to ficq — dD). In our 
second test case, this equation becomes: 

n(r,e») = npoic + rf!^(-^)%in''e (5) 

with R denoting the stellar radius. 

In this context, we can evaluate the total angular mo- 
mentum J within the convective zone by integrating over 
all elementary plasma torus of spherical coordinate (r, 6) 
and mass d?m (equal to 2tx p{r)r^ sm6 dOdr where p(r) is 
the local mass density, assumed to depend on r only). The 
resulting expression writes: 

»H /.7r/2 



J 



' sin 9 0.{r, 9) d m. 



(6) 



J 



dO) 



(7) 



Replacing S7(r, 9) with Eq. |1| and integrating over r and 
yields the following expression for J in our first test case: 

— mc R (iZcq - - 

where ttic represents the total mass in the convective zone 
and a the normalised mass-weighted second order moment 
of the fractional radius in the convective zone, given by 
rf (r/i?)^ dm/ dm (dm being the total mass of the ele- 

mentary plasma shell of radius r, equal to p{r)dr) . The 
same process applied to our second test case yields a similar 
expression for J, that writes: 

J mc R^ (ficq - A dQ.) (8) 

where 

A^l-i^ (9) 
5 a 

and (3 the fourth order moment of the fractional radius 
(r/i?)^ dm/ dm. These parameters mainly depend on 
the detail of the internal stellar structure and are listed in 
Table 01 for our three stars of interest, using the latest avail- 
able stellar structure models (Charbonnel et al. 1996; Siess 
et al. 2000). 

The conservation of angular momentum within the con- 
vective zone therefore implies, in our first test case, the fol- 
lowing relation between 57eq and dQ: 

1 



n. 



dO, + Qsb, 



(10) 



Q.ah being a constant, equal to the rotation rate the stellar 
convective zone would have if spinning as a solid body. In 
our second test case, this relation becomes: 



fie 



A dQ + Qsb , 



(11) 



implying in particular that the linear relationship between 
dn and fl^q is significantly less steep than in our first test 
case (see Table HJ- It may be worthwhile to note that, al- 
though f2eq and dfl refer specifically to the differential rota- 
tion parameters at the surface of the star. Eg. llOl still holds 
when replacing fl^q and dfi by their equivalents estimated 
at radial distance r from the centre of the star (and noted 
respectively f2eq(r) and df2(r)) since both do not show any 
dependence with r in this particular case; however, this is 
not true in the second test case, where Eq.^^b^comes, when 
expressed as a function of r2eq('') and dQ,{r): 



Table 4. Model parameters depending on the internal structure 
(sequentially the total mass M, the external radius _R, the frac- 
tional radius at the bottom of the convective zone rc/R, the frac- 
tional mass included in the convective zone rUc/M, the mean 
surface rotation rate (f2), the second, fourth and sixth order mo- 
ment of the fractional radius in the convective zone a, /3 and 7, 
as well as the model parameters A and fi defined in the text) for 
AB Dor, LQ Hya and HR 1099 



Parameter 


AB Dor 


LQ Hya 


HR 1099 


M (Mq) 


1.0 


0.95 


1.0 


R (Rq) 


1.05 


0.95 


3.7 


Vc/R 


0.73 


0.73 


0.16 


mc/M 


0.074 


0.074 


0.73 




1.4 X 10 


-* 4.5 X 10-5 


2.6 X 10-5 


a 


0.58 


0.58 


0.27 




0.35 


0.35 


0.11 


7 


0.22 


0.22 


0.054 


A 


0.52 


0.52 


0.68 


M 


0.027 


0.027 


0.035 


f2cq(r) = A'(r 


) dVL[r) 







(12) 



(13) 



where A'(r) is equal to 

with ^(r) denoting the squared fractional radius {r/R)"^. In 
our second test case, we therefore expect the relationship 
between d^{r) and r2cq(f) to change (turning steeper with 
decreasing r) depending on the depth at which the tracers 
used to estimate differential rotation are anchored. 

We can also derive by the same method an expression 
giving the kinetic energy K within the convective zone. The 
general formula, writing: 

/ 

Jtc Jo 

becomes, in our first test case: 



, 9)'^ d^m, 



(14) 



K — — mc 
3 



R' 



[Qcc - ^ dnf + ^ dfi' 



and in the second test case: 

K= - mc R^ [(ficq - A dQ,f + n dO?] 



(15) 



(16) 



where ^ stands for: 

M=-^-i^f^)' (17) 
^ 35 a 25 \a) ^ ' 

and 7 for the sixth order moment of the fractional radius 
\^ (r IRi)^ dm/ rf dm. In both cases, we obtain that K can 

JTc ^ic 

be written as the sum of two terms, one that we know is 
constant (given Eqs. 1101 and 1111 reflecting the conservation 
of angular momentum) and a second, potentially variable, 
one, directly proportional to dOF' . To compute the actual 
variation in kinetic energy dK, one must note that dO. can 
itself be written as the sum of two terms: 



dO. = dSln 



(18) 



where dfimin and cif2var respectively represent the minimum 
value of the observed photospheric shear (that correspond- 
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ing to the state closest to solid body rotation) and the vari- 
able contribution to dO,. In this context, the variation of 
kinetic energy induced by the differential rotation changes 
writes: 

dK = — rUc R fl (df2var + 2 d^var dQmin) (19) 

(/i being replaced by 8/175 in our first test case). 

It could be tempting to consider that the whole varia- 
tion in kinetic energy is transformed into magnetic energy 
as part of the magnetic cycle (as proposed by Lanza et al. 
1998). However, this conclusion first requires that we ob- 
tain an expression similar to that derived above for dK (see 
Eq. I19|l , but for the associated changes in potential energy 
that result from the fluctuations in the internal angular ve- 
locity field and thus to the variations in the shape of the 
star. We could then derive which fraction of the kinetic plus 
potential energy is converted into magnetic energy during 
the activity cycle. Similarly, we could also obtain estimates 
of the fluctuations in the in the stellar quadrupole moment 
(as done for instance by Applegate 1992), but in the specific 
context of our study (i.e. with the internal velocity fields 
assumed above in Eqs. |l]and|^. We however keep this out 
of this initial study and postpone it for future papers; it 
is nevertheless rather safe to claim that the amount of ki- 
netic energy transformed into magnetic energy is of the same 
magnitude as (though not strictly equal to) dK. 



4.3 Internal velocity fields and magnetic energy 
fluctuations of active stars 

If we now compare the model predictions obtained in 
Sect. 14.21 with our observations, we straightforwardly ob- 
tain that the second test case (featuring constant angular 
velocity along cylinders within the convective zone) seems 
much more appropriate. It is indeed only in this case that 
different depths in the convective zone correspond to dif- 
ferent angular shears (see Eqs. |1| and |^ , thus providing a 
natural explanation of the discrepancy between the flcq and 
dil estimates derived from Stokes / and V data. If this is 
indeed the correct interpretation, it would imply that mag- 
netic features, showing the strongest shear, are in general 
anchored closer to the surface (where dr2(r) is maximum) 
than cool spots. This is at least in qualitative agreement 
with the conclusions of, e.g., Donati et al. (2003), claiming 
that magnetic features (and in particular those hosting pre- 
dominantly azimuthal fields) must be formed very close to 
the surface in these stars. Further support about the idea 
that dynamo magnetic fields (and in particular their toroidal 
components) can be produced very close to the stellar sur- 
face comes from the recent discovery that the Sun features 
a subphotospheric layer with a strong radial gradient in an- 
gular rotation (Corbart & Thompson 2002) that seems at 
least able to trigger dynamo action (Dikpati et al. 2002). 

However, there is already a number of potential prob- 
lems with this simple description. It implies for instance that 
the range of possible latitudinal shears we should observe on 
a given star at a given epoch goes from dQ at the surface, 
down to dQ{rc) — dfl (,{rc) at the bottom of the convective 
zone, i.e. that the maximum ratio between the shears we 
measure from Stokes V and Stokes / data should be l/^{rc)- 
Although this is verified for AB Dor (for which the measured 



ratio varies from 1.25 to 1.8 and thus remains always smaller 
than the maximum predicted value of 1/0.73^ = 1.9), the sit- 
uation is different for LQ Hya where we observe (at epoch 
2001.99) differential rotation shears in Stokes V and Stokes 
/ that differ by a factor of about 6 (see Table|5]l, much larger 
than the theoretical maximum of 1.9. The origin of this dis- 
crepancy is not clear yet; the fact that only one point (out of 
the seven for which measurements from Stokes / and V data 
are available) exhibits this problem tells at least that the 
model we propose, although probably oversimplified, is not 
necessarily a complete nonsense. One possibility is that the 
model we consider does not take into account the additional 
shears that probably exist in the interface layer between the 
convective envelope and the radiative interior, nor that as- 
sociated with a potential subsurface shear similar to that of 
the Sun (Corbard & Thompson 2002); we may thus imagine 
that, in some rare cases, the bulk of the differential rotation 
tracers we use is anchored in these thin shear layers, pro- 
ducing a small number of observations in disagreement with 
our simple picture. It may also indicate that the differen- 
tial rotation estimate secured for LQ Hya at epoch 2001.99 
from Stokes I data is affected by some spurious contribution 
that we have not yet identified (values derived from Stokes 
V data being always more reliable and less prone to system- 
atic effects); although we are of course fully aware that is a 
rather unsatisfying and potentially worrying statement, we 
cannot yet rule out this possibility completely until a larger 
sample of estimates are available for us to check if this prob- 
lem is truly an exception or rather a regular occurence. In 
the rest of the study, we exclude this point from the analysis 
and further investigate how compatible is the toy model we 
propose with the rest of our data. 

A quantitative test consists in measuring the slope of 
the dfl versus 57oq relations that we obtain for each star 
by taking into account all estimates derived in Sect. |3 (and 
shown on Figs.|21to0J and comparing it to those the model 
predict (see Eqs. 1101 to 1121 . To determine the optimal lin- 
ear relationship that provides the best fit through all data 
points and their associated confidence ellipses for each star 
(something that a conventional least-squares fit, not taking 
into account the correlated errors on both axes, cannot prop- 
erly achieve) , we iteratively consider a large number of such 
relationships (i.e. a large number of A and ilsb pairs) and 
compute, for each of them, the closest approach (in terms of 
X^) to all error ellipses and the associated overall goodness 
of fit (defined as the sum of the minimum square distances 
to each ellipse). From these values, we derive a map (in 
the fisb — A plane), look for a minimum and fit the map 
by a bi-dimensional paraboloid in the neighbourhood of this 
minimum (by a method similar to that described in Sect.|5J 
to obtain both the optimum values of A and Q^h as well as 
the corresponding error bars. 

For AB Dor and using all the estimates listed in Tabled 
(10 points altogether with their associated error ellipses), 
we obtain from this procedure that the minimum achievable 
reduced is 10.5, already indicating that the model pro- 
vides a rather poor explanation of our observations. When 
2 (out of the 10 available) differential rotation estimates are 
removed from the full sample (those obtained from Stokes 
I data at epoch 1995.94 and from Stokes V data at epoch 
2001.99) on the reason that they deviate most from the fit 
including all points (at a level of about 5 a), the minimum 
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achievable reduced decreases to a smaller, but still un- 
comfortably high, value of 4.8, yielding values of A and Q,sh 
equal to 0.48 ± 0.08 and 12.212 ± 0.005 radd"^ respectively 
(with 1 G error bars enhanced by a factor of 2.2 = (4.8)'^''^ 
to reflect the poor fit obtained with this model). Note that 
the value of A we derive agrees well with that predicted for 
the second test model (equal to 0.52, see Table 0. If we 
now impose A to be equal to 0.2 (the value predicted for 
our first test model), the minimum reduced that we can 
now achieve is about 2.6 times larger than that obtained 
without constraining A. This strongly suggests that our first 
model (assuming solar-like differential rotation within the 
convective zone) is a much less probable option than the 
second one (assuming constant angular velocity along cylin- 
ders within the convective zone), event hough the accuracy 
to which the observations are fitted with this second model 
is still not satisfactory. 

To investigate whether the agreement with the data can 
be improved, the best way is to attempt exploiting all built- 
in characteristics of our simple model, for instance by tak- 
ing explicitly into account, when comparing with the data, 
the fact that the different types of tracers may be anchored 
at different convective depths. To achieve this, we proceed 
along the lines presented in Sect. 14.21 and try to measure 
the differential rotation estimates corresponding to Stokes 
I and Stokes V data with different linear relationships (us- 
ing Egs. 1121 and I13II . More specifically, we proceed as ex- 
plained in the paragraphs above, except that we now have 
three independent parameters, A (the slope of the f2cq versus 
dQ. relationship in the Stokes V data, presumably referring 
to the surface layers). A' (the slope of this relationship in 
the Stokes / data, presumably referring to deeper convec- 
tive layers) and one value of 0.sh- The only difference from 
the previous procedure is that we now fit a three dimen- 
sional (instead of a two dimensional) paraboloid to the 
map. This new modelling attempt (with the same two data 
points excluded) yields a minimum reduced of 3.2, along 
with values of A, A' and Slsb respectively equal to 0.34±0.06, 
0.28±0.08 and 12.222±0.004 radd"^ (error bars being again 
scaled up in proportion to the square root of the reduced x^)- 
Although the fit to the data is now more satisfactory and the 
derived value of A' falls within the expected range (between 
0.13 and 0.40, given the inverse squared fractional anchoring 
depths 1/5 obtained from ratioing the estimated shears at 
each epoch, see above), the value of A we extract is poorly 
compatible with that we expect for the surface layers (equal 
to 0.52), with a mismatch of about 3 a. Repeating the pro- 
cedure with A being fixed to the expected value of 0.5 yields 
no improvement in the accuracy level to which the data are 
fitted with respect to the two parameter fit described above. 

We therefore conclude that the estimates of differen- 
tial rotation and of its temporal variations at the surface of 
AB Dor suggest that the distribution of angular velocities 
within the convective zone is closer to that we expect for 
rapid rotators than to that of the Sun. However, the simple 
model we propose, in addition to the fact that the accuracy 
level to which it fits the data is only rough, also predicts 
some features that are not reproduced in the observations, 
and in particular the different slopes of the linear fieq versus 
dQ. relationships expected for tracers anchored at different 
convective depths. What it may simply indicate is that our 
test models are still far too simple to match the observations 



at that level of detail. Further attempts, with more realistic 
large scale angular velocity fields in the convective zone will 
be the subject of forthcoming papers. 

The first part of this experiment can also be applied to 
our LQ Hya data. Using the two parameter fitting proce- 
dure described above through the available measurements 
(4 points with associated error ellipses, once discarded the 
estimate from Stokes I data obtained at epoch 2001.99, see 
above), we obtain values of A and Slsb that are respectively 
equal to 0.36 ± 0.03 and 3.912 ± 0.003 radd"^ and provide 
a nice fit to the data, with an associated reduced x^ level 
of 0.84. For this star again, we obtain that the slope of the 
observed f2cq versus dQ linear relationship is not compatible 
with a solar like differential rotation within the convective 
zone; repeating the same procedure with A fixed to 0.2 (as 
required by our first test model) increases the minimum re- 
duced x^ up to 7.2, i.e. 8.4 times that achieved with the 
unconstrained fit. Assuming now constant angular velocity 
along cylinders, the model would then require the tracers to 
be anchored at about mid depth within the convective zone 
(at a fractional radius of about 0.86) to agree with the ob- 
servations. Another possibility is of course that the internal 
velocity field of LQ Hya is somehow intermediate between 
that of a very fast rotator like AB Dor and that of a slow ro- 
tator like the Sun, refiecting the fact that the rotation rate 
of LQ Hya (3.9 radd"'^) is bracketed by that of AB Dor 
(12.2 d) and that of the Sun (0.25 radd~^). In this context, 
contours of constant angular rotation within the convective 
zone should start deviating from being aligned with the ro- 
tation axis, and get slightly tilted towards the radial direc- 
tion. Once again, we caution that these conclusions are still 
highly speculative and require more observational material 
(in particular measurements on other stars rotating slower 
and faster than LQ Hya) to be settled more firmly. 

Carrying out a similar investigation in the particular 
case of HR 1099 would also be very interesting, as our second 
model predicts a significantly higher slope for the 57oq versus 
dQ. linear relation (equal to about 0.68, see Table [IJ due 
to the fact that the convective zone (pertaining more than 
99% of the total stellar volume) is much deeper in this star. 
The smaller rotation rate of HR 1099 (2.2 radd~^) may 
also act at the same time to decrease the value of A with 
respect to the predicted one. Since no temporal variations of 
differential rotation are however yet detected on HR 1099, all 
this remains essentially a matter of speculation, until more 
accurate estimates (such as those of Petit et al. 2002, but 
covering a longer time span) are available. 

Another important aspect of this problem concerns the 
amount of kinetic energy that the star either withdraws from 
or releases into the total energy reservoir of the convective 
zone when switching from one state of differential rotation 
to another. In particular, we should make sure that the to- 
tal power associated with the energy transfer from or into 
kinetic energy is significantly smaller, or at the very least 
not larger, than the stellar luminosity itself, to ensure that 
the process we invoke to explain our observations (which in 
principle should apply to the vast majority of rapidly rotat- 
ing cool stars) does not implicitly require at the same time 
that the star looses its energy at a very fast rate. Of course, 
we reckon that the energy transfers we invoke likely cor- 
respond to some redistribution process between the various 
reservoirs in which the convective zone can store energy (i.e. 
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kinetic, potential and magnetic energy wells); however, we 
suspect that these transfers are inevitably associated with 
some energy dissipation, i.e. that a significant fraction (the 
exact proportion being poorly known) of the energy trans- 
fered is lost (and thus radiated) in the process. 

The changes in kinetic energy associated with the dif- 
ferential rotation fiuctuations can easily be estimated from 
Eq. [IHl In the case of AB Dor, the results reported in 
Sect. l3.l1 suggest that dfimin and dfivar are respectively equal 
to 60 and 40 mradd"^, implying that dK should be equal 
to Aafi nicR^dflvar^ /3 or about 4 x 10^^ erg. Assuming this 
transfer occurs progressively over a typical timespan of a few 
years (10* s), we obtain that the power required to operate 
the change we invoke corresponds to about 1% of the total 
stellar luminosity. (Note that dK is only weakly dependent 
on the assumed internal velocity field, varying only by a fac- 
tor of about two between our two test models). For LQ Hya, 
our observations indicate that dSlmin is probably close to 
zero while dSlvar gets as large as 200 mradd"^, implying 
that dK should now be equal to afj, mcJ?^df2var^/3 or about 
2 X 10*° erg. If we again consider that the corresponding 
evolution occurs on a time scale of a few years, we find that 
the associated power corresponds to 10% of the total stellar 
luminosity. Although this sounds rather large already, one 
must keep in mind that only a small fraction of this amount 
is actually lost (i.e. radiated) in the process; it suggests that, 
even in the case of the star that exhibits the most drastic ef- 
fect (i.e. LQ Hya) , the corresponding power dissipated in the 
process does probably not exceed 1% of the stellar luminos- 
ity. Since no temporal changes were detected in the case of 
HR 1099, we thus cannot obtain a similar estimate without 
any (hopefully reasonable) assumption on dilmin and dSlvar- 
Taking dilmin = and df^var = 40 mradd"^ as a guess (and 
to remain grossly compatible with the differential rotation 
changes that the orbital period fluctuations observed on this 
system request, Donati 1999), we obtain that dK reaches as 
much as 7 X 10*" erg and that the associated power required 
to drive this change corresponds to about 3% of the stel- 
lar luminosity (provided the change in dQ occurs over a few 
years) . 

If we finally assume that most of the changes in kinetic 
energy associated with the fluctuations in differential rota- 
tion are associated with equal (and opposite) variations in 
magnetic energy occurring throughout the whole convective 
zone, we can obtain a very rough estimate of the average 
magnetic field strength required to drive these changes. We 
caution that this estimate is necessarily very rough for at 
least several reasons. The first reason is of course that the 
magnetic field is expected to be very intermittent through- 
out the convective zone and to exhibit a strong gradient 
in the radial direction (comparable to the one that the gas 
pressure undergoes), implying that an average magnetic field 
strength over the whole convective volume is not extremely 
meaningful; the second reason (already mentioned above) is 
that fluctuations in potential energy are also expected to oc- 
cur (as a consequence of the changes in the shape of the star) 
and probably participate as well in the total energy budget. 
However, a rough estimate being better than no estimate, we 
decided to mention in the paper the derived average mag- 
netic strengths, to at least convince the readers that the 
magnitude of these fields is physically plausible. For AB Dor 
and LQ Hya, we obtain magnetic intensities of 9 and 23 kG 



respectively; using the above guesses in the case of HR 1099 
yields a field strength of 5 kG. This is clearly larger that 
the magnetic strengths estimated at the surface (e.g. Do- 
nati 1999; Donati et al. 2003), and reflects the strong radial 
gradient of field strength (pointing inwards) that probably 
exists within the convective zone. 



5 CONCLUSION AND PROSPECTIVES 

This paper reports positive detections of surface differential 
rotation on two rapidly rotating cool stars at several epochs, 
by using stellar surface features (both cool spots and mag- 
netic regions) as tracers of the large scale latitudinal shear 
that distorts the convective envelope in this type of stars. 
We also report definite evidence that this differential rota- 
tion is different when estimated from cool spots or magnetic 
regions, and that it undergoes temporal fluctuations of po- 
tentially large amplitude on a time scale of a few years. 

We consider these results as further evidence that the 
dynamo processes operating in these stars are distributed 
throughout the convective zone rather than being confined 
at its base as in the Sun. With the help of two very simple 
models of the angular velocity field within the convective 
zone (one that resembles that of the Sun and one supposed 
to represent that of ultra fast rotators), we explore the quan- 
titative implications that angular momentum conservation 
within the convective zone imposes on potential variations 
of surface differential rotation with time. 

This modelling suggests, from comparison with our ob- 
servations, that the rotation velocity field within the convec- 
tive zone of the stars we investigated is not like that of the 
Sun, and may resemble that we expect for rapid rotators, 
with constant angular velocity along cylinders aligned with 
the rotation axis. The power required to drive the observed 
changes in the differential rotation of the convective zone are 
estimated to be of the order of 1 to 10% of the total stellar 
luminosity; we speculate that these changes result from the 
dynamo processes (and from the underlying magnetic cy- 
cle) that periodically converts magnetic energy into kinetic 
energy and vice versa. 

The technique outlined in this paper corresponds to the 
first practical method for investigating the large scale rota- 
tion velocity field in convective zones of cool active stars. In 
particular, we note that it can already provide us with more 
information than what one could obtain with other exist- 
ing techniques such as asteroseismology (expected to yield 
rather poor resolution on the internal rotation profile in the 
radial direction, even when operated from space, at least in 
the next decade). We therefore strongly argue in favour of 
renewing the experiment reported here, in order to explore 
more extensively the full capabilities of this new method. 
Monitoring continuously the same star for time spans as 
long as several months every year for several years would 
be especially interesting if one wants to disclose short term 
chaotic variability (presumably due to the natural random- 
ness inherent to turbulence and differential rotation) from 
the long term supposedly periodic evolution driven by the 
magnetic cycle. 

In this context, the new spectropolarimeter ESPaDOnS, 
designed for maximum efficiency and being built for the 
3.6 m Canada-France-Hawaii Telescope (commissioning 
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scheduled for spring 2003), represents the optimal instru- 
mental facility for carrying the observations that this new 
method requires. We expect in particular that the large gain 
in sensitivity that this new tool will provide over older gen- 
eration instruments should give us the opportunity of both 
extending the available data set to a much larger sample 
of stars and studying the physical processes of interest to 
us in much more details. This research program will also 
strongly benefit from having similar instruments on smaller 
telescopes to allow the long, repeated and if possible mul- 
tisite, observational campaigns that are necessary for opti- 
mising the performances of the investigation we suggest. In 
this respect, the 2 m Bernard Lyot Telescope atop Pic du 
Midi in France, that will soon be equipped with NARVAL (a 
carbon copy of ESPaDOnS, whose construction just started) 
and dedicated to spectropolarimetric studies, should be ex- 
tremely useful for this kind of research. 



Petit P., Donati J.-F., Wade G.A., ct al., 2003b, MNRAS (sub- 
mitted) 

Press W.H., Teukolsky S.A., Vetterling W.T., Flannery B.P., 
1992, in: Numerical Recipes in C, Cambridge University 
Press, p. 689 

Reiners A., Schmitt J.H.M.M., 2002, A&A 384, 155 
Siess L., Dufour E., Forestini M., 2000, A&A 358, 593 
Strassmeier K.G., Bartus J., 2000, A&A 354, 537 
Vogt S.S., Hatzes A.P., Misch A.A., Kurster M., 1999, ApJS 121, 
547 

Zahn J.-R, 1989, A&A 220, 112 



ACKNOWLEDGEMENTS 

We warmly thank Corinne Charbonnel for kindly providing 
us with some models of the internal structure of the three 
stars on which this study is focused. We are also grateful 
to the referee, Rob Jefferies, for making very detailed and 
relevant comments that allowed us to improve the paper 
substantially. 



REFERENCES 

Applegate J.H., 1992, ApJ 385, 621 

Brown S.F., Donati J.-F., Rees D.E., Semel M., 1991, A&A 250, 

463 

Brunning D.H., 1981, ApJ 248, 274 

Charbonnel C, Meynet G., Maeder A., Schaerer D., 1996, A&AS 
115, 339 

Cameron A.C., Donati J.-F., 2002, MNRAS 329, L23 
Cameron A.C., Donati J.-F., Semel M., 2002, MNRAS 330, 699 
Corbard T., Thompson M., 2002, Solar Phys. 205, 211 
Dikpati M., Corbard T., Thompson M.J., Gilman P.A., 2002, ApJ 
575, L41 

Donati J.-F., 1999, MNRAS 302, 457 
Donati J.-F., Brown S.F., 1997, A&A 326, 1135 
Donati J.-F., Brown S.F., Semel M., Rees D.E., et al., 1992, A&A 
265, 682 

Donati J.-F., Cameron A.C., 1997, MNRAS 291, 1 

Donati J.-F., Cameron A.C., Hussain G.A.J. , Semel M., 1999, 

MNRAS 302, 437 
Donati J.-F., Cameron A.C., Semel M., ct al., 2003, MNRAS 

(submitted) 

Donati J.-F., Mcngcl M., Carter B.D., Cameron A.C., Wichmann 

R., 2000, MNRAS 316, 699 
Hall D.S., 1990, in: Ibanoglu C. (ed.), "Active Close Binaries". 

Kluwer, Dordrecht, p. 95 
Howe R., Christensen-Dalsgaard J., Hill F., et al., 2000, Science 

287, 2456 

Donati J.-F., Semel M., Carter B.D., Rees D.E., Cameron A.C., 

1997, MNRAS 291, 658 
Johns- KruU CM., Valcnti J.A., 1996, ApJ 459, L95 
Kitchatinov L.L., Riidiger G., 1999, A&A 344, 911 
Lanza A.F., Rodono M., Rosner R., 1998, MNRAS 296, 893 
Petit P., Donati J.-F., Cameron A.C., 2002, MNRAS 334, 374 
Petit P., Donati J.-F., Wade G.A., et al., 2003a, MNRAS (sub- 
mitted) 



